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1. Introduction 


Analysis of electrically large structures 





Q Multiple excitation analysis of electrically large structure 
= Geometry to be analyzed fixed 


" Scattering characteristics analyzed according to 
electromagnetic source change. 


> Required for analysis of reflection characteristics in each 
direction of structures 





<Example of monostatic RCS for multiple direction> 





Rigorous methods Asymptotic methods 


EE = Rigorous method rather than asymptotic method 
" Strength in scattering analysis 


= Only boundary of problem is meshed 


Advanced Computational ; 
Electromagnetic Laboratory http://www.ausairpower.net/APA-2011-03.html 





(GI) YONSEI UNIVERSITY 


1. Introduction 


Method of Mom ent 


Ani 


Integral equation > Matrix equation > Matrix solver 
Z l= V, Z : Geometry characteristic matrix, | : Unknown current, V : 


Electromagnetic source 


O(N?) memory, O(N?) CPU time required, N : Number of unknowns 


Two types of matrix solver 
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e Strength in memory/CPU time for single 
excitation 
« ZI z in : Current distribution for n, 
iteration 
* Independent computation for different V 
* Time-consuming for multi-excitation problems 
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Reduction of DoFs using macro basis functions 
Free from convergence problems 

Strength in multiple excitation problems 
Single Z-t computation for different excitation 
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1. Introduction 


Characteristic basis function method 





Q CBFM 
" Applicable for arbitrary geometry 


= Good compression effectiveness (Number of MBFs / Number of low-level basis functions) 


k 


Pa " Create macro basis 





rn —À functions 
a (Characteristic basis) 
uu. -— à that reflects the 
- — . e 
Sa characteristics of each 
1 NS block 
" Divide geometry into blocks eN 


Making impedance matrix for CBs |: = Solution of reduced matrix equation 


o : " Reorganize current distribution based on CB to 
zw =< J zo Jio > current distribution based on low-level basis 
b1b2 =< Jb1; 451 42) 7 b2 





aracteristic BasissFunction Method. In: Computational Electromagnetics. Springer, New York, NY, 2014. p. 1-40 
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1. Introduction 


Acceleration of CBFM 





" Making impedance matrixfor CBS | a Most time-consuming procedure in CBFM 





VR =< Joi. ze Jio > > Acceleration scheme using multipole expansion (MCBFM) 


" Trade-off in MCBFM | =" Applying multipole expansion to 
| multilevel CBFM (MLCBFM) 
= When block size is small ; = Various block size for multiple level 
v Sampling rate : " Proposal of a technique with merits 
Y Nestronerato : for each block size 








= When block size is large 


Y” Compression rate 





v Number of transfer functions : 
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2. CBFM and MCBFM formulation 












Excite planewaves 
Number of planewaves relates 
with bandwidth of radiation 
spectrum of block 

Now > 2(kr,+275)* 

20x20x2 (Theta x phi x TM/TE) 
for cube edge length 0.5 2 





2 V 


= Divide geometry with specific block size | 7 

TRUNK MU : mxn mxm mxn nxn 
= Extended region for fictitious edge effect : 
= 0.5 2 À (Cube edge length) ~ Singular value decomposition process 


* 0.14 (Extended region) | = Elimination of redundant CBs which have singular value 
Advanced Computational | lower than threshold 
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2. CBFM and MCBFM formulation 


CBFM formulation 





" Reaction terms of EFIE (RWG basis) 





= Excitation terms of EFIE 


" Reaction terms between CBs 





Ny No = Representation of CB 
zi = Fu: LFranddS = f ER - LD 1 ia as Na 
S a ; | Fim = ` DE 
Noi Noo ) (n) A NEP p 
— ` Du. In Phi ` L(fp2,5) dS > E =< Ji, ae Jp > 


= Excitation terms for CBs 





VO = | Fi, (r) - E/^*(r E insu - E"*(v)dr 
> jo J f(x) - E/"*(r > Jy > FWO 
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2. CBFM and MCBFM formulation 


CBFM formulation 





" Reduced matrix equation 


eG S cA Rs ee ehh [LE] [ein 
E A < h 2; hA» e «A20. || EP < Ja ums 
< Jp, 2. Jy > < Jp, 2) gre ves BINA SS n | c da VO s 





; SVD : 0.1 h (0.4%) 
Other : 1.1 h (3.6%) CB gen mat fili : 0.4 h (1.396 
lm. ^^ gen mat inv : 0.4 h (1.396) 


Reudced mat inv : 7.8 h (25.496) Composition 


Ratio of 
Sean Time MoM mat gen : 8.5 h (27.396) 


Analyzing Aircraft 
@ 400 MHz 
with 170,000 Unknowns 






Reudced mat gen : 12.6 h (40.5%) 


= Bottleneck process of CBFM 
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2. CBFM and MCBFM formulation 


MCBFM formulation 





" Multipole expansion technique for Green's function 


Maii jkke( Ty (k, k, ra )dS 
AME" NA Ea Wed LE mod 
tlr — r'| | idi 
= nd k) Ty (k. k Tap) T(r’, k)dS 
1 
where 





= 


ij 
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2. CBFM and MCBFM formulation 


MCBFM formulation 





= Reduced matrix filling with addition theorem 


e iklr-r| 
Yay IO [jwp J fins: J fj» ;(1 - pV) der) 
No1 Nba . 
M D J R(f;, k) - Tr (k, k, ra) T(f;, k)dS 
where 


R(f;, k) = jun | f;(r)e -jFkerra qy. 


2 


O-level Aggregated/Disaggregated CB 





T(f;, k) — [1 — kk]. Je re) Pero dy! 
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2. CBFM and MCBFM formulation 





MCBFM formulation 


= Reduced matrix filling with addition theorem 


Ny Noo 
um J (XO T? R(f;. k)) -Tr (k, k, ra) (O 155), T(E;.k))dS 
i j 


_ J R(Fotm:k) - Tr (k. k, rab) T (Foon, k)dS 
1 


where 
. Ny | 
RAF him: k) = y Di RE, k) 
No 1-level Aggregated/Disaggregated CB 
T(Fozn, k) 2 V^ I] T(E, k) 
j 


= Complexity of reduced matrix filling : O(N¿pp*Np?) 
= Complexity of reduced matrix filling using multipole expansion: O(Ncpr?) 
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3. Multilevel MCBFM 


Multilevel CBFM 





= A sequence of blocks concept is defined 


= [seq] = [i,,i,,...,i,], where the levels go from 
the highest to lower ones. 


" (Characteristic basis for high level 
Biseq] Niseg],b 


E [seg),p - 2. y D b, t dim 


TIL 


" (Characteristic basis for lowest level 


Nseg] 
E seg], m 2 Jo [seq], 7 s 


= Reduced impedance matrix for level / 





(1) 2T 0-1) = 
sca; [seg] ^ € " Iseai]^ A [seg] [segj]: “ [seas] > 
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3. Multilevel MCBFM 


P | PP. 
a E * \/ | À - 
| ; = L uw a uU \ A Y À 
CHO ADA AD A CX D | E À 4 Jii m ff 
"— 
- 





= Three types of interaction zone 


" Near-zone interaction 
Y Near-zone in high-level block 
v Near-zone in low-level block 





= Intermediate-zone interaction 
Y Near-zone in high-level block 
V Far-zone in low-level block 


= Far-zone interaction 
Y Far-zone in high-level block 
V Far-zone in low-level block 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Near-zone interaction 


(2) | . 
LA) [Bra > / Era] p * Ein] g)d5 


Nia}, Nip 
=f (p) | | 
o J >, Irati mE m 29» a „FIB n)dS 

ME Nana 


| (p) (q) 7 
Ri LO las Finds 


Ny A],2 Nig; i 


T / >. lom Faam ` L( X IBF un)dS 
Ni A].2 Na 
du 2. Hp FA dom  L (> I9 a ,Fip.215)d$ 
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CBs in each block 


Nia] A 


Nia? 


Fr = DE 1: mF(4,1]m +1 om FA aj 


NÑNiB],1 


NiBj2 


| (p) (p) 
FiB] p = » m 1mY[B,1],m +2 La. > mE (BY), a 


Tr 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Intermediate-zone interaction 


Nipy,1 


¡DAA (Ta 19 F1B1)m)d5 


Nia Nia 


= À 2, n im E bbs | Fa} * L(F(B1jn)dS 


Nia Niet - 
x = 2, 2, I8. Lum Ds Bam | Raam 
Tr (kk ria] — rea) T (Fan k)ds 


NA], 


= [i >. LE à RCE [A 1, k)] 


Ng. 
Tr (E k raa] — riga V I9 1 T (Fia) K)]dS 
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CBs in each block 


Frajp = 


Fiap 


ÑNta],1 Nij, 


DE naim +1 2m F IA 2m 


Niet T NiBj2 ió 
= > lim Bm + 2 amina 


TTL 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Intermediate-zone interaction 


= [Rí (Frasi Kk) Tr (E k rp — rie) (Fina: k)dS 
where 
[A],1 
R(F - > F 14,1}, m: k) 
Nia | 
T(FiB] 1n K) B > I ın TFB k) 
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CBs in each block 
Niaj Niaj.2 


DE naim +1 2m F IA 2m 


Frajp = 


Nin, - N181,2 a 
Fiap = >, Ip Lm FB tje +2 Li 2. T [B,2], Tn 


TTL 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Far-zone interaction 


(2) 
Z AJ.[C],pa — J, Frajp L(Ficjq)d5 


Niaj, 1 E 1 ( 
- 2. n "nem m'L i Clin Ficin)ds 


TL 


NiA)4 2 


Nie}, 
+] 2. BA Lm F [A, 1],m ' A, IE nas 


CBs in each block 








Nia], Ni a],2 
N Nin 3 
P IV af = I ] o E [4] p = PO 2, I 1. mE A, 1] |,m + >. i E mi (A, 2],m 
> yl 2, mi A,2],m ( >. | iC] 1,5 E [C11,5) 05 m 
Nic ,1 Ntc], J 
NiA],2 Ntc1,2 e o -Y I) | IC m » m Eo Fc]. v 
+ 2, lá IA am F Aa LA 3 T1652, F1215)d5 m 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Far-zone interaction 


Ny A],1 Nic] ,1 


B 2. > Idea J, Fini: Lied 


Tr TL 


Niaj, Nic],2 


NI o . 
F Y >, Hum lien | Fin» | L(Fic21n)d5 


Nj41,2 Nic],1 


Er (p) (a) 
T 2. > Hamlin | Fram: L(Frcjn)ds 


Nraj,2 Nie, 


(n r(q) ge | . 
* 2. y T amie 2n | F ra 9m  L(Fio21,5,)d5 
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CBs in each block 
Niaj, Nia? 


TN P) 7 (p) 
FiA] p = ` I 1 E [A]m + >. bhi 2.m¥ [A,2],m 


Tr 


Nicj,2 


| (D) 4 (p) 
Fic]p = >, I6, 1 s FIca]m + >, L101,2mE[C,2,m 


TIL 


Ntc] 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Far-zone interaction 


Nia. a Nies 


"2 > 14 ],1,m Prenant 


Tr (E, k, ria] — ry) T (Fio 


—jkk-(riA 1]—r[A]) 
ke gkk- (rc. i] rie) ds 


Niaj 1 Nic] 2 


| (p ikk-(r;, —r[ 
«Y X fth. ef [oe om De Aa a 


TL(k, kr — ra) T(Fcan; K)e jk (rica -rie) qs 


Nia].2 Nic] „1 


+ 2 2 T am Ii | RE ay - ke —jkk-(ri4 j—r(4]) 


TL, kr — rg) T(Fiea s ket (rica ren qs 


NL A], 2 Nic].2 
(q) | ns db ri E 
* » L I VERLO J R(FyA 3, kje aaa) 


TL (k, k, ria] — rpg) (Fia. 
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k)e jkk- (ric 21 ric) ds 


= CBs in each block 
Fia]p — Po Y yi 1. mE A, 1],m + y UA E mt DE 2],m 


Nie},1 Nic},2 


== >, LG} mF calm + » TO 2m € [C2] m 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Far-zone interaction 


= Ji Y Tijm RC Atm kk e 3E eria 17a] - Ty (k, f, r4] — r[c]) 
Nie}, 1 
> D ln T(Fic 1]. Kk] i^e (Tic, 1] "el S 
Ni]. a | 
+ | >, D dm CE [A 1]. a k)e _jkk- (Tha, 1] T [A] E -Tr (k, k, r LA] T ric) 
Nias 
| 2 ña T(Fic, 2].n Ke jkk (Tie, 2] rici ]dss 
NA], 2 


+ ff b If 7 2m RCE {A 3, Uu k)e e IUE -TL(k, k,  T[A] — r[c]) 


X ne ln T(Fic,1],n: k)e jkk: (Tic,1] reds 


MEE 
+ fl Le 17 2m RCE A 2m k)e 7 tnt) - Tr (E, K, ru] — rie) 


Nie ].2 


2 Le, ln TFC in: 
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= CBs in each block 
Nia}, Ni A)2 


Fia] p = 2, IFA. m t 2, i 2 2 mE [42m 


Ny Ct Ni C1,2 


| (P) 4 (p) 
CLp >, LG} mF calm + » TO 2m € [C2] m 
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3. Multilevel MCBFM 


Multilevel MCBFM 





= Far-zone interaction 


Nia 


p / [DS Lam RCE [Ass Ke” Peta rap 
Naja —— 
+ D a Rai De Pt) 
Tr (k, k, F[A] — ric) 
Nicis 


DETTES (Fic, i, K)e/ e rie rie) 


Nie}.2 


E D Hea TE cam inetrearnemas 


- / R(F a) pk) Tr (E, K ria] — rg) T(Fqoj 4, K)4S 
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CBs in each block 


Fiap = 


Njaj,1 Ni A)2 
Do AmE t D 2m" [A,2],m 


Nic], Nic],2 


| (PD) (p) 
E >, LG} mF calm + » TO 2m € [C2] m 








m , BS 
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4. Numerical results 


Analysis of 3 A radius sphere 





Q Proposed 


Mie series 


| | | | | | M- 2-level 
Time [s] CBFM 
120,000 | | 


, CBFs generation 1,134 
es 3,558 


© 
D A wf 


RCS [dBsm] 


DR PE OÙ D 
EU X AD A 
= S) ) 

S 9 (D 


RR RR TRY rp en AA AD 
VV: n Y: € + 6 

AN NA NN SEC S6 

SAR LES RS RR 

a C 

QU 


Reduced matrix 7.351 
filling 





0 20 40 60 80 100 120 140 160 180 


0 [degree] 
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5. Summary 


Summary 





OW) Algorithm for multiple excitation analysis of electrically large 
structure is proposed 
" Proposal of a technique with merits for each block size 
= Based on CBFM, MCBFM and MLCBFM, algorithm has been proposed 


Q Three types of interaction calculation scheme and multilevel 
aggregated/disaggregated CBs 
= Near-zone interactions, intermediate-zone interactions, far-zone 
interactions are computed for different way 


" High level aggregated/disaggregated CBs are derived through 
reorganization of low level aggregated/disaggregated CBs 


Q Numerical results shows good agreement and efficient 
computation of electrically large structure 


Advanced Computational 


Electromagnetic Laboratory e 





(Ed) YONSEI UNIVERSITY 


Thank You for Your Attention 


Advanced Computational A 
Electromagnetic Laboratory 


